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Abstract 

We report a study of disordered electron systems with spin-orbit coupling on a cylinder using 
methods of random matrix ensembles. With a threading flux turned on, the single particle levels 
will generally avoid, rather than cross, each other. Our numerical study of the level-avoiding gaps 
in the disordered Rashba model demonstrates that the normalized gap distribution is of a universal 
form, independent of the random strength and the system size. For small gaps it exhibits a linear 
behavior, while for large gaps it decays exponentially. A framework based on matrix mechanical 
models is suggested, and is verified to reproduce the universal linear behavior at small gaps. Thus 
we propose to use the distribution of the spectral-flow gaps associated with flux insertion as a new 
way to characterize 2d random systems with spin-orbit coupling. The relevance and qualitative 
implications for spin (Hall) transport are also addressed. 
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Introduction In recent years, there is a surge of interests in doing spintronics, i.e. to ma- 
nipulate and to control electron spins in solid state devices, with external electric, instead 
of magnetic, fields. Potential advantages include efficient generation of spin polarization 
inside the devices, handy manipulation of individual spins at nano-scales, and desired re- 
duction of heat dissipation during information processing. (The last point is because, in 
particular, the electric field can induce a transverse spin Hall current, respecting time rever- 
sal invariance and therefore being dissipationless Since there is no fundamental direct 
coupling between spin and the electric field, the spin-orbit coupling (SOC) plays a central 
role in inducing the spin transport. For these reasons, the study of spin transport and, in 
particular, the spin Hall effect in systems with SOC has recently attracted much attention 
(for recent reviews, see e.g. 2|, yj). In the literature, a mechanism for the spin transport due 
to SOC is called extrinsic or intrinsic, depending on whether disorder (such as impurities 
or imperfections) is involved or not. However, in a real sample there is always disorder; 
to study possible interplay between the intrinsic and extrinsic mechanisms, it is better to 
study them simultaneously in one unified framework, rather than to study them separately 
and add their contributions. In this paper we report such a study of the systems with SOC 
using the methods of random matrix ensembles. 

A prototype of the intrinsic SOC is the Rashba coupling in a two dimensional electron 
gas. The intrinsic contribution [4] to the spin Hall conductivity in this model, which has a 
universal value c/Stt, was found to be exactly cancelled by the extrinsic vertex contribution 
from impurities in a perturbative approach 0, ^. Naturally one asks whether the cancella- 
tion persists beyond perturbation theory, and if so whether this points to a new universality 
to be discovered in the disordered systems with SOC concerning the spin (Hall) transport. 
These are the questions we want to address. 

We put the Rashba model with on-site randomness on a cylinder, and mimic an infinites- 
imal electric field by turning on slowly a magnetic flux threading the cylinder. The latter is 
equivalent to imposing a twisted phase in the boundary condition. The same technique has 
been used before in studying charge transport in the presence of localization [3| , as well as 
the integer Sj and fractional [9] quantum (charge) Hall effect. What we are concerned about 
here is the behavior of level-avoiding gaps due to disorder as the threading flux varies. (The 
relevance of the level-avoiding gaps to the transport has been discussed before in Refs. 1^ 



for the charge Hall effect and in Refs. UJ, ll3|] for the spin Hall effect in this and similar 
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models.) We have carried out a numerical study and indeed found a new universality in 
the random system with SOC. More concretely, the distribution of the (normalized) level- 
avoiding gaps is found to have a universal form, independent of both the random strength 
and the system size. It is linear at small gaps and decays exponentially at large gaps; and 
the average gap for un-normalized distribution seems to tend to a finite value when the 
system size becomes very large. In this paper we will suggest a matrix mechanical model for 
the underlying random ensemble, which is verified to exhibit the desired small-gap behavior. 
Usually one uses the level-spacing statistics as a characteristic of the universality in a ran- 
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15| on a square 



dom system or ensemble Here our results suggest a different type of characterization 
of random ensembles using the gap statistics in the spectral flow due to flux threading (or 
twisted boundary conditions). 

The Disordered Rashba Model Consider the disordered Rashba model 
lattice with size L x L: 

H = -tJ2 4cj + ^iAci + V'so U XI - i Y cjayd+s: + h.c. J , (1) 

(ij) i \ i i / 

where q = (ci|,Cj|) are electron annihilation operators for up-spin and down-spin, respec- 
tively, t the nearest-neighbor hopping, and Vso the Rashba spin-orbit coupling. Disorder 
is represented by uniformly distributed on-site random potential Wi G [—W/2,W/2]. We 
impose the free boundary conditions in x-direction and the twisted boundary conditions in 
y-direction: \E'(i, j-|-L) = e*^\E'(i, j), so that the model is put on a cylinder, which is threaded 
by a magnetic flux ^ = h6/e. 

We diagonalize the Hamiltonian numerically. The basic parameters are L, W and 
Vso- We have taken their ranges to be L=12, 14, 16, 18, 20, 24, 28, 32, 36, 40, 48 and 
W = 0.05t,0.1t,0.2t,0.4t,0.6t,0.8t, l.Ot, with the SOC fixed to be Vso = 0.5t unless stated 
otherwise. We can safely assume that our system is in the metallic phase, as the range of the 
randomness, W, is much less than the critical value (Wc = 5.875t for Vso = 0.577t) fl5^ for 
the insulating phase. We have first calculated the level spacing distribution, Phsi^',^) with 
6 fixed, where s is the spacing between adjacent levels. Indeed we have seen a flux driven 
crossover from the symplectic to the unitary ensemble. Some typical results are shown in 
Fig. [T]for = and 6 = 1. We find the following power law behavior in the small-s region, 

, s^ for 9 = 0, symplectic, , , 

PLs{s;e)<x{ (2) 
s for 6 = 1, unitary. 
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FIG. 1: The level spacing distributions phs{s;0) for = and 9 = 1. L = 24 and W = 0.2t. For 
both data, the numbers of the level spacings we take are about 10^. 

a) b) 




FIG. 2: Level diagrams for L = 10 a) without disorder (W = 0) and b) with disorder (W = O.lt). 

in agreement with the random matrix ensembles with the corresponding symmetry jl^ : 
When 6 = 0, the system has time-reversal symmetry, leading to Kramer's degeneracy, which 
is broken when 6' 7^ 0. 

Spectral- Flow Gap Distribution We are mainly interested in the spectral flow due to 
a varying flux. Typical level-flux diagrams are shown in Fig. [2l Without disorder (i.e. 
W = 0), the energy levels cross each other as 9 is varied (Fig. [2^), because momentum ky 
is a good quantum number. With W 0, however, the translation symmetry in y-direction 
is broken by disorder, and generically gaps open near the would-be level-crossings (Fig. 
Wp)- To characterize the spectral flow, one might try to extract from numerical data the 
statistical distribution of the level-avoiding gap g which is, however, expected to depend on 
the parameters of the model. So a better way to detect whether there is universality is to 
plot the distribution pciO of the normalized gap ^, the value of g normalized by the average 
{9)'- i = 9/ {9) with L and W flxed. 
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FIG. 3: a) The gap distributions Pg(?) foi' various L and W (in units of t). The soUd curve is the 
function ([3]). b) The gap distribution Pg(0 for L = 24 and W = 0.2t. For comparison, the function 
([3]) (solid hne) and the Gaussian unitary ensemble (GUE) surmise {iT^/2)e~'^^^^^ {igI (dashed line) 
are also shown. The inset shows the small-^ behavior which can be fitted by a linear behavior oc ^ 
(dashed line). The number of gaps counted for each data point is greater than 5 x 10^. 

Our numerical data for a wide range of L and W, as exemplified in Fig. [3^, demonstrates 
a remarkable universal behavior, namely the distribution, pciO^ of normalized gaps 
neither depends on L nor W. Indeed the distribution is fit nicely by the following function: 

= ^i^' ""^kWy wwy 

where ({z) is the zeta function, and constants a and b are fixed by the normalization condi- 
tions /q°° dC^PciO = 1 /q°° d^^pQ^^) = 1. In Figl3)D, we show the data for the behavior of 
the distribution at large and small gaps: It is linear at small gaps and exponentially decays 
at large gaps. 

To examine whether the SOC is crucial for the appearance of the universal behavior or 



not, we present in Fig H] some data for the variance, v = y d^^'^pciOi of normaUzed 
gap distribution with and without the SOC, respectively. While the data with Vso = 0.5t 
and Vso = 0.25t support a universal behavior, those with Vso = clearly show non-universal 
L- and ly-dependence. 

Numerically we have also checked the large-L behavior of the average gap {g), which is 
known to be non-universal. Figure [5^ shows that it tends to a non-zero value as L ^ oo 
with W fixed. This fact is essential for the distribution of the normalized gap ([3]) to make 
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FIG. 4: The variances Va^ of the gap distribution a) without and b, c) with SOC. For comparison 
the variance of the function ([3])(sohd hne) and that of the GUE surmise (dashed hne) are also 
shown. 
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FIG. 5: a) Non-universal behaviors of the average gaps {g). b) The VF-dependence of the average 
gap in the limit L — > oo. 

sense in the thermodynamic limit. When combined with the Laughhn's argument, it has a 
profound imphcation concerning spin Hall transport, which we are going to discuss below. In 
FigEb, we show the I^-dependence of the limiting value {g)o; the solid curve is the following 
function, (^)o = at ta.nh {(3\n{-fW/t + 1)) with a = 2007.8, /3 = 0.50 and 7 = 7.50. 

New Universality of Random Ensembles The above-discovered universality of the distri- 
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bution Pg{0 inspires us to search for the random ensemble underlying it. We notice that 
the spectral flow gaps we focus here are similar to the so-called avoided crossings that have 
been extensively studied in quantum chaos (see, e.g. l?!). Our model contains a periodic 
(flux) parameter 6 in the twisted boundary condition. Since time reversal invariance is re- 
spected only for = and vr, our system should correspond to a random ensemble with 
unitary symmetry except at ^ = and vr, near which there are crossovers to symplectic 
symmetry. The distribution of avoided crossings for Gaussian unitary ensembles (GUE) was 
studied before in the literature of quantum chaos 17||. In Figl3)D, we show the comparison 
of our data with the GUE surmise given in Ref. IGj. Though the linear behavior in the 
small-^ region agrees with the GUE studied in Refs. [l^, the exponential decay in the 



large-^ region certainly does not 



16| . Since the spectral flow is due to a varying parameter 



20l |. in which the 



9, a more natural setting should be a matrix quantum mechanical model 
time variable t may be identified with 6. Below let us show that in a wide class of matrix 
mechanical models, the distribution pg(0 is indeed linear at small gaps and is independent 
of either the size of matrices or the potential in the model. The Lagrangian of the models 
we consider is 

L = -Tr$(t)2 + TiVmt)) , (4) 

where $(t) is an x Hermitian matrix, and TrV^($) a potential invariant under the 
automorphism $ G^G^ with G unitary. Writing $(t) = f/(t)Miag(Ai(t), ■ • ■ AAr(t))f/(t) in 
terms of its eigenvalues Aj(t) and the unitary matrix U {t) that diagonalizes it, the Lagrangian 
is recast into 

^ = E + E 1^^ - I' ^^^^n + E ^(^^) (5) 

i i<j i 

with uj = iUU"^ . The corresponding Hamiltonian is 

^ = ^Ep.^ + Em) + En;^^, (6) 

^ i i i<j 1^* " 

where pi and Lij are the conjugate momenta of Aj and ujji, respectively. The statistical 
properties are determined by the Gibbs measure, 

dGN = ^e~^^ n dkdp, n dLi.dL*^, (7) 
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FIG. 6: An avoiding gap between Ai and A2. 



where Zj^ is the normahzation constant. The un-normahzed gap distribution n(^g) is given 
by 

N 



{g) = iV(iV - 1) n / dG^ - tn)S{-g/2 - \,)5{g/2 - A2), (8) 

fc=3"'!?/2<|Afe| „ 



where tn denotes a level avoiding point between Ai and A2. See FigjGl From the normalization 
condition, the gap distribution is obtained by 

(9) 



Pg(0 = {9)n{,g) / dgn{g) 



At the level avoiding points t = tn, we have pi = p2 with P2 — Pi > 0. Thus J2n ^(t ~ ^n) in 
the integral can be replaced with 



Y,Kt-tn) = 5(P2 -Pl)(P2 -Pl + |P2 -Pl|)/2 



-S{p2-Pl] 



OH 



OH 



diX2 - Ai) 



(10) 



5(A2-Ai) 

where we have used the equations of motion for pi and p2- While the resulting integral is 
hard to perform, the smaAl-g behavior can be easily seen to be 



nig) = Ag, 



(11) 



where A is a constant given by 
4iV(A^-l) ^ 



A 



x5{p2 - Pi) exp 



N N N 

II ^dp^l[ [ du, [ dL*^ n /rfAfc n (^-^ - ^^)' n ^ 

i=l '' i<j •' k=3 2<k<l k=3 

N 



N 2 

I E I - E - E ^(^fc) - 21^(0) 

i=l ^ j<j fc=3 



(12) 



Thus we have a universal hnear behavior for pciO the small-^ region for any N >2. The 
insensibihty of the hnear behavior to the number of levels coincides with that to L for 
the Rashba model as the number of levels of the Rashba model is given by 2L^. 

We note that the large-.^ behavior of the random matrix quantum mechanics does depend 
on the form of the potential V{X). This is well-known for usual random matrix ensembles. 
For example, for polynomial potentials, the kernel of the unitary ensembles in the large N 



limit has a universal form corresponding to the Gaussian ensemble 



2ll: 



sinTrfx — -u) , . 

K{x,y =^ f^; 13 

7r{x - y) 

but for logarithmically growing potentials, V{\) ~ ^(InA)^ (A ^ 1), it shows another 



universality 



a sinvrfx — y) ,^ 
vr smh a\x — y) 



It is plausible that there exists a class of potentials in matrix mechanical models, which 
reproduce the spectral-fiow gap statistics we have found from our numerical calculations. 

Spectral- Flow Gaps and the Hall Transport In the above we have studied the distribution 
of spectral-flow gaps in a disordered system with SOC, as a response to twisted boundary 
conditions. This is in line with the Laughlin gauge argument jsl for the charge Hall effect, 
which has been recently adapted to the case of the spin Hall effect [ll, 12, 13 1 for systems 



with SOC. Our data show no level crossing in the spectral flow, in agreement with Refs. 
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13| for the disordered models. Because of level avoiding, each level always flows back to 
its original energy after the flux adiabatically increases from zero to the unit flux quantum. In 
contrast to the quantum (charge) Hall effect, in the present case the system is in the metallic 
phase, and there are levels going across the Fermi level. At absolute zero temperature and 
in the adiabatic limit, the levels that are below the Fermi level and never go across it will 
maintain its own contribution to the spin Hall conductivity, similar to the case of the charge 
Hall effect. The sum of these levels to the spin Hall conductivity of the system is expected 
generically to have zero expectation value if the position of the Fermi level is set randomly. 
However, during the adiabatic variation of the flux, the instantaneous spin Hall conductivity 
will fluctuate because of the levels that move out of or dive into the Fermi surface. So at 
least for the cylindrical geometry, the time-averaged spin Hall conductivity should be very 
small, compared to the clean intrinsic limit. 



Moreover, a finite electric field will induce the transition across the Fermi surface between 
levels with sufficiently small gaps (Landau-Zener-like tunneling). Also at finite temperature, 
thermal fluctuations will cause similar transitions between levels with gaps smaller than 
the average thermal energy. Therefore, the fluctuations in spin Hall conductivity may be 
signiflcant if the temperature is not too low and the applied electric fleld is not too small. 
It is no doubt that the knowledge we obtained in this paper on the distribution of level- 
avoiding gaps and on the magnitude of the average gap will be useful for detailed analyses 
of the disorder effects on the fluctuations in spin Hall conductivity. Here we only point 
out some features inferred from our results. First, the level avoiding gap distribution is of 
a universal form for the normalized gap, which is linear at small gaps. So we expect that 
there is also a universality in the fluctuations in a^"^ if properly scaled. In other words, the 
fluctuations may have a one-parameter scaling law by the average gap at sufficiently low 
temperature, where the Landau-Zener tunneling is dominating. Second, since the average 
gap is non- vanishing in the large-L limit, we expect that the fluctuations in o"^™ will be 
governed by an exponential factor involving the temperature or the electric fleld over the 
average gap. 

Summary and Discussions In this paper, we examined the distributions of level avoiding 
gaps in the Rashba model in the presence of disorder. We found that the normalized dis- 
tribution is of a universal form, independent of the random strength and the system size, 
while the average value of the gaps is non- universal. 

A natural expectation is that our universal distribution works for other models with 
SOC. (For example, a disordered two-dimensional electron gas with Dresselhaus SOC [231], a 



genera, 
gas 



ized Rashba model 



24| . a graphene model with SOC 25|, and a two dimensional hole 



26|.) We note that this does not mean the spin Hall conductivity itself is universal: It 
depends on the average value of the level-avoiding gaps, thus it could be model dependent. 
For example, if the average value becomes zero in the thermodynamic limit, a net spin Hall 
effect can be observed. In addition, what we discuss in the spectral flow argument is not 
the spin Hall current itself but the spin accumulation after the unit flux insertion. Because 
the spin is not conserved in general, these two quantities are not necessarily in accordance 
with each other. This also supports the expectation that systems with different spin Hall 
conductivities may have the universal gap distribution we discovered. 

In the absence of disorder the Rashba model has integrability of the momentum con- 
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servation which is the origin of the level crossings. The integrability may explain a mixed 
behavior of our distribution between the Poisson distribution and the GUE one. Indeed, it 
is known that the gap distributions for fully chaotic systems without integrability show no 
mixed behavior and are in good agreement with the GUE surmise 

To conclude, we remark that the new universality we discovered in this paper may have 
consequences concerning quantum transport in mesoscopic systems with SOC. It may also 
provide helpful insight into the interplay between localization and SOC in the spin Hall 
effect of disordered systems. 
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